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The physical phenomenon of interest is high-speed gas dynamics 
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Aerodynamics 
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Properly characterizing the aerodynamic performance of reentry 
vehicles is critical for optimal trajectory design. 
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heat transfer rate — » thermal protection material selection 
heat load — > thermal protection material thickness 



Governing Equations 
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per unit mass, and P is the pressure. 



The viscous stress tensor r and the heat flux vector q are defined as 
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Equation (4) with (7) is known as Stokes’ hypothesis for a Newtonian 
fluid [2], 



Provided the transport properties may be related to the 
unknowns, Equations (l)-(3) form 5 equations in 7 unknowns 



is enlarged to contain 

a local mass balance statements for constituent species 
o local energy balance/exchange statements for internal modes 



In the literature equations (l)-(3) are often treated as the system 1 
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the notation is cumbersome, but it is fairly standard 
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Upwinding is required to stabilize convection-dominated flows. In SUPG 
schemes this is accomplished by biasing the test functions W in the upstream 
direction: 
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For compressible flows t supg is generally diagonal [3]. 



SUPG stabilization does not yield monotone solutions. Additional treatment 
is needed to prevent spurious oscillations in regions of shockwaves. 
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The shock capturing operator, S, was adapted for a system of conservation 
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In this work the numerator is modified from its original form to retain 
consistency with (9) in the case of transient, viscous flows. 



Expand U(x, t) and F,(x, t ) in terms of standard Lagrange finite element basis functions: 
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Expand U(x, t) and F,(x, t ) in terms of standard Lagrange finite element basis functions: 
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Expand U(x, t) and F,(x, t ) in terms of standard Lagrange finite element basis functions: 
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contrasts to the typical approach in which 



Consider a one-dimensional, inviscid, normal shock at Mach 5. For this 
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which implies that pu, pir + P, and puH are all constant. 



Inviscid Flux Discretization 
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The semidiscrete weak form in equation (13) is discretized in time using a backwards finite 
difference scheme. Both first and second-order accurate in time schemes may be derived from 
Taylor series expansions in time about U„+i : 
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Linearization 
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Newton Scheme 




Newton Scheme 
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Newton Scheme - Comments 
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This approach is especially attractive for “real gas” flows 



The resulting sparse implicit linear systems are amenable to solution with 



This may be approximated via differencing residual evaluations, yielding a 
so-called matrix-free method: 



Implicit Solution Strategies 
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This work uses matrix & matrix-free 
GMRES with full ILU-0 preconditioning 
linearization is important 
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The X- 1 5 Experience 




Edney’s Type IV Interaction Pattern [10] 



An experimental test program was conducted in 1998 by France’s Office National d‘Etudes et 
de Recherches Aerospatiales (ONERA) to investigate shock-shock interactions produced by an 




Static temperature contours 
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Type IV Shock Interaction 



0 (degrees) 



Hypersonic flow over a missile nose tip with a forward facing cavity has been observed to 
exhibit transient flowfield response in both experimental investigations and numerical 
simulations [12, 13]. The flowfield response characteristics are largely driven by the cavity 
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investigate the presence of 
vibrational nonequilibrium in the 
freestream [161 
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Steady states, runs 2893 and 2894 




Time Convergence, run 2894 



Time Step 



High speed schlieren, run 2890 





Computed schlieren, run 2890 



Possible Mechanism for Observed Unsteadiness 
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Noise Characterization [17] 
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Variation of Pitot Pressure Fluctuation 
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Results - Flowfield 





Results - Surface Pressure 



Normalized Distance From Cone Apex, 



25 kHz, 6% RMS Pitot Pressure Fluctuation 



Nondimensional Distance From Cone Apex, 



Frequency Influence 



Normalized Distance From Cone Apex, x/D 
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